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Signatures of minor mergers in Milky Way-like disc 
kinematics: Ringing revisited 
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ABSTRACT 

By means of A^-body simulations we study the response of a galactic disc to a 
minor merger event. We find that non-self-gravitating, spiral-like features are induced 
in the thick disc. As we have shown in a previous work, this "ringing" also leaves an 
imprint in velocity space (the u-v plane) in small spatial regions, such as the solar 
neighbourhood. As the disc relaxes after the event, clumps in the u-v plane get closer 
with time, allowing us to estimate the time of impact. In addition to confirming the 
possibility of this diagnostic, here we show that in a more realistic scenario, the in-fall 
trajectory of the perturber gives rise to an azimuthal dependence of the structure in 
phase-space. We also find that the space defined by the energy and angular momentum 
of stars is a better choice than velocity space, as clumps remain visible even in large 
local volumes. This makes their observational detection much easier since one need 
not be restricted to a small spatial volume. We show that information about the time 
of impact, the mass of the perturber, and its trajectory is stored in the kinematics of 
disc stars. 

Key words: galaxies: formation - galaxies: kinematics and dynamics - methods: 
analytical - methods: A^-body simulations 



1 INTRODUCTION 

The present day spatial distribution, motion, ages 
and chemical abundances of stars contain impor- 
tant information about how gal axies form and evolve 
IIFreeman &: Bland- HawthornI [20021 ). Substructure in these 



or self-gravitating 
[Fux '2001'; 'Minchev 



spiral arm s (see e.g. Dehnen 2000l : 
Quillcn '2008"; jAntoia et al] l2009l : 
Minch ev ct al. 2010; Quillcn ct al. 20111)" or debris from 



Way (e.g. Ibata. Gilmore & Irwin 1994: Helmi et al. 


19991; 


Newberg et al."200^; 'ibata et al. 2003; Yanny et al] 


20031; 
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2009: 


Helmi et al.l 20061; Starkenburg et all 120091; iKlement et all 


20091; Williams et all |2009|. 201 ll), but also in 


other 


galaxies fe.E. Ibata et al.l 200ll; Chapman et al.l 


20081; 


Martmez-Dekado et al.l |2008|. |2009|; iMcConnachie et all 


20091) . This substructure can arise from a variety of different 



host fe.E. iHelmi & White 
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20081; Michel-Dansac et al 


201 ih. Another Dossible mecha- 


nism. introduced by Minchev et al.l l|2009l. hereafter M09). 



physical mechanisms. Examples are the disruption of open 
clusters (|Eggenl 1 19961 ). the dynamical effects induced by 
non-axysimmetric components of a galaxy, such as a bar 
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is related to the response of a galactic disc to the tidal 
interaction with a satellite galaxy as it merges. Thanks to 
dynamical friction, relatively massive satellites can reach 
the centre of the disc before being completely disrupted 
by the tidal field ( Binney fc Tremaine 2008 ). M09 showed 
that the sudden energy kick imparted by the gravitational 
potential of the satellite, as it crosses the plane of the disc, 
can strongly perturb the velocity field of disc stars located 
in local volumes such as the Solar Neighbourhood. These 
perturbations can be observed in the u-v plane as arc-like 
features travelling in the direction of positive v, where u and 
V are the radial and tangential velocity components. This 
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mechanism, known as "ringing" , could explain the presence 
of high velocity streams observed in the Solar Neighbour - 
hood, such as the Arcturus s tream I Williams et al.|[ 2009h . 
Furthermore, as shown by iQuillen e t al, (2009), mmor 
mergers can also induce radial mixing in the outer disc, 
thus partially explaining the lack of correlation between 
age and metallicity observed in the Milky Way's thick disc 
fsee I Villalobod [20091 . and references therein). 

An important limitation in the M09 analysis is the ex- 
tent of the local volume that can be explored to search for 
signatures of ringing in the u-v plane. As the depth of the 
volume is increased, a progressively larger number of den- 
sity waves is included. This causes the waves to interfere, 
wiping out the substructure previously observed in this pro- 
jection of phase-space. This not only places constraints on 
the size of stellar samples that can be explored to search 
for signatures of ringing in the Solar Neighbourhood, but 
also makes studying this mechanism in fully self-consistent 
A'"-body simulations unfeasible. Note that to resolve local 
volumes as small as, e.g., 100 pc, an extremely large parti- 
cle resolution is required. Thus, in M09 only relatively sim- 
ple high-resolution 2D simulations of massless test-particles 
were performed to explore this problem. 

In this paper we show that the space of energy and an- 
gular momentum is much more suitable for identifying the 
signatures of ringing in a galactic disc. In this projection of 
phase-space much larger local volumes can be explored with- 
out erasing the substructure associated with ringing (Sec- 
tion [3]). This allows us to analyse a set of fully self-consistent 
A'^-body simulations of a merger between a relatively mas- 
sive satellite galaxy and a pre-existing thin disc (Section Q. 
We briefly describe our simulations in Section[2l In Section[5] 
we summarise our results and discuss its application to cur- 
rently available as well as forthcoming full 6D phase-space 
samples of thick disc stars in the Solar vicinity. 



2 SIMULATIONS 

Two types of numerical simulations are studied in this work. 
First, we analyse a two dimensional massless test-particle 
simulation of an unrelaxed stellar galactic disc, as described 
by M09. In this model, a fiat rotational velocity profile is 
considered, with a value of 220 km s~^ at 7.8 kpc from the 
galactic centre. The disc follows an exponential density pro- 
file with a scale-length of 3 kpc. Initial velocities for the par- 
ticles are drawn from Gaussian distributions in the u and v 
directions. In particular, the standard deviation in the radial 
direction is set at au = 50 km s"'^. To emulate a population 
that is unrelaxed, velocities are purposely chosen so that 
particles are unevenly distributed in their radial oscillation. 
A total of 7 X 10^ test-particles were used to represent the 
disc. For more details, we refer the reader to M09. 

Second, we analyse a set of fully self-consistent N- 
body simulations of the merger between a satellite galaxy 
and a pre-existing galactic thi n disc. Most of these sim- 
ulations were first presented by IVillalobos fc Helmil l|2008l . 
hereafter, VH08) to study the formation of a galactic thick 
disc due to the heating induced by the merger event. 
In particular, the simulations considered in this work are 
those referred as "z = 1" in VH08, which were set up 
to emulate a merger event that could have given rise to 



Table 1. Initial properties of the three A'^-body simulations analysed 
in this work. 



Host's DM halo 


Virial mass, Mf^l 

' host 

Virial radius 
Concentration 


5.07 X 10" 

122.22 

6.56 






[A/©] 
[kpc] 


Host's stellar disc 


Mass, Aff 
Scale length, i?d 
Scale height, zq 


1.2 X 10^° 

1.65 

0.165 






[Me] 

[kpc] 
[kpc] 


Satellite's DM halo 


/ = 20% 


15% 


10% 




Virial mass, M^^f^ 
Virial radius 
Concentration 


1.01 

71.35 

8.09 


0.75 

64.65 

8.41 


0.5 

56.63 

8.86 


lO"[A/0] 
[kpc] 


Satellite's stellar bulge 


Mass, A/=^Jr"="' 
Scale radius, rb 


2.4 
0.709 


1.8 
0.654 


1.2 
0.583 


W^[Mq] 
[kpc] 



the Milky Way's thick disc. The models for the host and 
satellite galaxies consist of two fully self-consistent compo- 
nents: a dark matter (DM) ha lo following a NFW profile 
(jNavarro. Frenk fc Whit3ll996i '). a stellar disc for the host, 
and a bulge for the satellite system. The stellar disc fol- 
lows an exponential density profile, and the satellite's bulge 
a Hernquist profile (|Hernauistlll99Gl ). The properties of the 
Milky Way-like host and the sate llite systems were scaled t o 
those expected at z = 1 follow i ng IMo. Mao fc Whit^ (|l998l ). 
As described bv lBrooks et al.l (|201l[ ). in this picture gas in a 
halo conserves its specific angular momentum (equal to that 
of the dark matter) as it cools to form a centrifugally sup- 
ported disc that grows from the inside out. In the simplest 
model, in which the density profile of the galaxy is modelled 
as a singular isothermal sphere, the radius of the resulting 
disc is proportional to the parent halo virial radius. Thus, 
in a A-CDM cosmology, this would imply a decrease in the 
disc size of a factor of 1.7 out to z = 1. Later studies have 
shown that the evolution of the disc size - stellar mass re- 
lation predicted by this model is considera bly stronger than 
what observations are suggesting (see e.g. ISomerville et al.] 
'2008!). Note however that the results presented in this work 
are not strongly dependent on the particular value chosen 
for the host disc's scale-length. 

Both disc and bulge are represented with 1 x 10^ A''- 
body particles. To study the impact that mergers of different 
total mass ratios have on the phase-space structure of the 
disc, three different values of / = Mi';°st/Afsat were explored. 
These are / = 20%, 15%, and 10% . Table [T] summarises the 
values of the parameters involved in these models. Note that 
/ also sets M^^'^JM^^ and Mllfl^^/Mf^''^^ 

The following is a brief summary of the evolution of the 
merger simulations analysed in this study. We refer the in- 
terested reader to VH08 for a complete description. In all 
cases the satellite was released from a distance of 84 kpc 
(« 50 -Rd) from the galactic centre in the most likely in- 
falling orbit according to previous studies on distr ibutions 
of or bital parameters of infalling substructure (e.g. iBensonI 
l2005h . Here Rd is the host disc's scale-length. In such rather 
eccentric orbit (e = 0.86), the satellite rapidly reaches the 
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central region of the galactic halo where it is more effec- 
tively affected by dynamical friction and tidal disruption. 
This causes the satellite to lose orbital angular momentum, 
which progressively shrinks its infalling orbit. In such a way, 
for the simulation with / = 20%, 15%, and 10%, the satel- 
lite takes roughly tmerger El 2, 2.5 and 4 Gyr to reach and 
settle in the disk centre, respectively. In particular, for the 
/ = 20%, the satellite experiences ~ 8 pericentric passages, 
all of them at less than 20 kpc from the disk centre. As the 
satellite decays in its orbits, it induces the formation of tidal 
arms in the disc, which carry angular momentum from the 
disk centre outwards. Once the satellite is brought on to the 
plane of the disc by dynamical friction, it transfers kinetic 
energy from its orbits to the disc stars, increasing their ver- 
tical motions and causing a visible thickening. By the end 
of the simulations the morphological, structural and kine- 
matical properties of the thickened disc and the disrupted 
satellite settle down and become stable, which takes place 
roughly 2 Gyr after the satellites have reached the disk cen- 
tre. The final scale-heights of the resulting thick disks are 
« 5 zo for / = 20% and « 4 zo for / = 10% (see Figure 14 
of VII08). Here zo is the scale-height of the pre-existing thin 
disc. 



3 TEST-PARTICLE SIMULATIONS OF AN 
UNRELAXED DISC 

In this section we briefly review the main results presented 
in M09 and study how the observed substructure in the ve- 
locity space of local volumes is distributed in E-Lz space. 
In the top panels of Figure[T]we show, at six different times, 
the distribution in the u-v plane of particles located within a 
ring of 160 pc width centred at 7.8 kpc from the disc centre. 
Note that the simulation is axysimmetric. Thus, the rings 
can be regarded as local volumes. As a consequence of the 
sudden energy kick imparted on the initially relaxed disc, a 
radial oscillation is excited on the disc particles. The radial 
frequency of these orbits will depend on both their orbital 
energies and angular momenta (see M09). Thus, particles 
with different energies will have different radial frequencies. 
These differences in frequencies induce phase-wrapping in 
velocity space: at a given time after the perturbation, in a 
given volume of space, we will observe groups of particles 
that have completed a different number of radial oscilla- 
tions. Each one of these groups will have a characteristic 
frequency or energy. These groups can be observed in the 
u-v plane as distinct arc-like features. Note that, as time 
passes by and the disk relaxes, the arcs get closer together 
and their number increases. The bottom panels of Fig. [T] 
show the same distribution of particles now projected onto 
E-Lz space. Particles that previously were part of a given 
arc in the u-v plane are now distributed along well-defined 
features of nearly constant energy. Due to the wide range 
covered in the v velocity component, each arc is stretched 
along the Lz direction. Note that, within a volume as small 

^ In this work the merger time, tmcrgcr, is defined as the time 
when the satellite runs out of orbital angular momentum. For the 
kind of orbits considered in this work, the oscillations between 
the centres of mass of both systems ceases approximately at the 
same time. 



as the one considered here, the gravitational potential, "I>(r), 
can be regarded as a constant. Thus, for groups of particles 
with approximately constant energy, 

= 2{E - $(r)) f» est, 

which explains the arc-like shape of the features observed in 
the u-v plane. 

As pointed out by M09, due to the differential rotation 
of the galactic disc, the separation between arcs in the u-v 
plane depends on the galactocentric distance of the volume 
under study. Increasing the size of this volume, and thus 
sampling a larger range of galactocentric radii, causes the 
waves to interfere, wiping out most of the observed sub- 
structure in velocity space. This places strong constraints 
on the size of the sample of stars that can be explored to 
identify signatures of ringing in the solar vicinity if we use 
velocity space. 

In Figure [2] we compare the effect of volume depth in 
the u-v and E-Lz planes for the same simulation and time 
step. The top row we show the velocity distribution of par- 
ticles located inside rings with different widths. The mean 
galactocentric distance rmcan of each ring is 7.8 kpc and 
the maximum distance sampled from rmean by these rings 
is specified in each panel as dmax. Note that, by including 
in our sample particles that are only at 250 pc away from 
7'mcan, the wavcs Start to interfere and the arcs in u-v to 
overlap. In a ring with dmax = 2.5 kpc all the substructure 
has been completely erased. The bottom panel of Figure [2] 
shows the distribution of the same sets of particles in E- 
Lz space. When very small local volumes are considered, 
i.e. dmax = 0.08 kpc, both spaces are equivalent in terms 
of the observed substructure. However, it is very interesting 
to see that even in the ring with largest dmax = 2.5 kpc 
substructure can be clearly observed in E-Lz space. By in- 
cluding particles at larger radii we include density waves at 
different energy levels. In addition, each wave is better de- 
fined since particles that were in a different phase of their 
oscillation, and thus outside the arc, are now included in 
the sample. Note that, as expected, the separation between 
density waves decreases as a function of energy, i.e., with 
galactocentric radius. 

The results presented in this section are important for 
two reasons. First, by sampling relatively large local vol- 
umes, it is possible to study the properties of ringing in 
fully self-consistent Ai'-body simulations, where the parti- 
cle resolution is much lower than in the previous exper- 
iment. Second, future astrometric missions such as Gaia 
(|Perrvman et al.ll200ll ) will provide very accurate measure- 
ments of 6D phase-space coordinates for stars located at 
distance s as large as 2.5 k pc from the position of the Sun 
(see e.g. iGomez et al]|2010t ). Although the u-v plane and the 
E-Lz space are equivalent when searching for signatures of 
ringing in local spheres with radii as small as 80 pc, the frac- 
tion of thick disc stars in this volume will be considerably 
smaller than in a sphere of 2.5 kpc radius. Note that it is in 
the thick disc where we would expect to find signatures of 
ringing, since these stars are less affected by perturbations 
from the Galactic bar, spiral structure or giant molecular 
clouds. Thus, the E-L^^ space will allow us to explore much 
larger samples of thick disc stars from, e.g., a future Gaia 
catalog. 
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Figure 1. Distribution of particles located within a ring of 160 pc width, centred at 7.8 kpc from the galactic centre, obtained from our 
test-particles simulation of an unrolaxed disc. The different colours (contours) indicate different number of particles. The first column 
of panels shows the distribution of particles before the "merger" event. Top panels: Time development of the u - v plane distribution of 
particles. Note that, as explained in M09, as time passes by the number of arc-like features increases and get closer together. Bottom 
panels: Time development of the E-L^ distribution of particles. The substructure present in the u - v can be now observed as group of 
particles with approximately constant energy covering a wide range of angular momentum. 
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Figure 2. Distribution of particles within rings centred at 7.8 kpc 
from the galactic centre, obtained from our test-particles simula- 
tion of an unrelaxed disc. The different colours (contours) indicate 
different number of particles. Each column is associated with a 
ring of a different width, as indicated in the top panels. All dis- 
tribution are drawn from the same snapshot, t = 0.69 Gyr. As 
the width of the ring is increased, the arcs in the u-v plane (top 
panels) become less clear. At dmax = 2.5 kpc substructure has 
been completely wiped out. Note, however, that in E-Lz space 
substructure can be clearly observed. As expected, the number of 
features in this space increases with increasing dmax- Note as well 
that the separation between energy waves decreases as a function 
of energy (see text). 

4 FULLY SELF-CONSISTENT iV-BODY 
SIMULATIONS 

We now explore fully self-consistent TV-body simulations 
of minor merger events. As described by VH08, the spa- 



tial and kinematical structures of the pre-existing thin 
disc are strongly affected by the merger event. As the 
satellite decays in its orbit, it passes through the plane 
of the disc several times. At each passage, disc parti- 
cles are impulsively heated by the satelli te's potential 
(see e.g., iQuinn. Hernauist fc Fullagaij [1993 . VH08). One 
manifestations of this heating mechanism is the forma- 
tion of transient, non-self - gravitating spir al arms (see e.g. 
'Tutu kov fc Fedorova|[2cio^ : ioh et al.li2008h . The right panel 
of Figure [3] shows a contour plot of the distribution of par- 
ticles in the X - Y plane of our / = 20% simulation, 1 Gyr 
after fmcrgcr- It is very interesting to compare this distribu- 
tion with the one obtained from the test-particles simula- 
tion, shown in the left panel of the same figure. Since the 
energy kick imparted on the test-particles is axisymmetric 
annular features (left panel) rather than spiral arms (right 
panel), can be observed. It is important to emphasise that 
the spiral features observed in our A^-body simulations are 
not self-gravitating, as they quickly disappear with time. 
This can also be inferred from the large values of the radial 
velocity dispersion of the resulting thick disc shown in Fig- 
ure 3] In this figure we show the radial velocity dispersion 
profile of the disc at three different times. The radial veloc- 
ity dispersion is considerably larger at all radii already 0.1 
Gyr after the satellite galaxy is fully merged. This profile 
does not evolve significantly during the next 2.9 Gyr of sim- 
ulation. Within 5 to 10 kpc from the disc centre, the radial 
velocity dispersion takes a value of (Ju ~ 50 km , which 
is consistent with the value used in the test-particles sim- 
ulation. As expected, the asymmetric drift induced in this 
way causes a decrease in the rotational velocity as a function 
of radius after the merger. Note that, although its magni- 
tude decreased by about 20% at i? ^ 4 kpc, a flat profile 
is preserved. We will now analyse solar-neighbourhood-like 
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Figure 3. Left panel: Distribution of particles in the X - Y 
plane, as obtained from our test-particles simulation of an un- 
relaxed disc, ^ 1 Gyr after the energy kick. The different colours 
(contours) indicate different number of particle. Note the multi- 
ple annular features, associated with density waves. Right panel: 
Distribution of particles in the X - Y plane 1 Gyr after tmerger (t 
= 3 Gyr), as obtained from our A'^-body simulation of a merger 
event with a mass ratio of / = 20%. As a consequence of the non- 
axysimmetrical energy kick imparted by the satellite as it crosses 
the plane of the disc, spiral features (rather than annular) are 
formed. 



volumes in the resulting thick disc, to search for signatures 
of ringing in energy and angular momentum space. 




R [kpc] 



Figure 4. Radial velocity dispersion and mean azimuthal veloc- 
ity profiles at three different times, from our / = 20% Af-body 
simulation. The profiles have been computed within concentric 
rings of 1 kpc width for particles with |Z| < 1.5 kpc. The solid 
and dashed lines show the profiles computed at t = 2.1 and 4.5 
Gyr, which corresponds to 0.1 and 2.5 Gyr after 

^merger , respec- 
tively. The dotted dashed lines represent both profiles at t = 
Gyr. Note that the later were obtained after evolving the host for 
3 Gyr in isolation. The red vertical lines indicate the positions of 
the spheres studied in Section |4. II 



4.1 Density waves in E-L^ space 

Since the particle resolution of these simulations is much 
lower than that of our test-particles simulation, it is not fea- 
sible to look for signatures of density waves in volumes as 
small as the ones analysed in Section |31 In what follows, 
we focus our analysis on samples of particles located within 
spheres of 2.5 kpc radius. The simulations analysed in this 
work were designed to emulate the result of a minor merger 
experienced by a Milky Way-like galaxy at z = 1. Therefore, 
the scales of our final thick disc are relatively smaller com- 
pared to those obse r ved in the Milky Way. As explained by 
IVillalobos fc Helmil (|2009l '). it is not straightforward to de- 
cide which radii corresponds to a solar neighbourhood-like 
volume. Following their work, we place our sphere at two 
different locations, namely 5.5 kpc ('-^ 2.4_Rtd) and 8 kpc 
(~ 3.5i?td), where iitd = 2.28 kpc is the final scale length 
of the simulated thick disc. This also allows us to study the 
expected dependance of the number and distribution of den- 
sity waves with galactocentric radius. 

The left panel of Figure [S] shows the distribution in the 
u — v plane of particles that are located inside a sphere cen- 
tred at (8,0,0) kpc from the galactic centre, 0.7 Gyr after 
tmergeJEI. This distribution was obtained from our / = 20% 
simulation. Note that only 865 disc particles are found 
within this sphere. Due to the large volume sampled, the 
absence of arc-like features in this projection of phase- 
space is not surprising. The same distribution of parti- 
cles projected in E-L^ space is shown on the right panel 
of this figure. Interestingly, particles are now distributed 
in well-defined clumps of approximately constant energy. 

^ This particular time and location were chosen arbitrarily for 
simplicity. Dependencies with azimuthal angles are explored in 
Section 112] 



Each clump represents a different density wave crossing the 
sphere. Figure [6] shows the time evolution of this distribu- 
tion. Note that our solar neighbohood-like volume is rotat- 
ing with an angular frequency set by the velocity of the 
local standard of rest (LSR), obtained after correcting the 
final mean rotational velocity profile for axysirnmetr ic drift 
l|Noordermeer. Merrifield fc Arag6n-Salamancall2008l ) . Over 
time, the number of clumps increases and they come closer 
together, as expected from Section |3] (see also M09). Sub- 
structure can be clearly observed already at t « 0.8 Gyr and 
can be followed until t ~ 4.7 Gyr, spanning a total of about 
4 Gyr. Later, although still present, the signatures of the 
density waves are difficult to observe even in E-L^ space. 
Recall that for this simulation, tmerger ~ 2 Gyr. Thus, sub- 
structure can be observe in the E-Lz plane for about 2.7 Gyr 
after the satellite has completely "sunk" to the centre of the 
host's disc. Figure [7] shows, as a function of time, the E-L^ 
distribution of the particles located inside the sphere centred 
at 5.47 kpc from the galactic centre. Again, a wealth of sub- 
structure can be observed. Comparison of this figure with 
Figure [S] reveals the presence of a larger number of clumps 
at any given time after tmorgor- This refiects the dependance 
of the phase-mixing time-scale on galactocentric distance 
(see M09). Note that in this inner sphere the perturbation 
is triggered after ~ 1.1 Gyr, a time that corresponds to the 
second pericentric passage (and thus, closer to the galactic 
centre) of the satellite galaxy. Substructure can be clearly 
observed for a total of ~ 2.8 Gyr, and for ~ 1.9 Gyr after 

^merger ■ 

We have also analysed a self-consistent simulation of a 
galactic disk evolved in isolation. The disc is first relaxed 
within a rigid version of the halo potential for about 1 Gyr. 
Once the disc is relaxed the rigid halo is replaced by its 
N-body counterpart and the system is evolved fully self- 
consistently for 4 Gyr (See VH08, Section A3 for a detailed 
description of this procedure). Even after 4 Gyr of evolu- 
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tion, no noticeable substructure can be observed in the E-Lz 
distribution of particles located within solar neighbourhood- 
like volumes. 



4.2 Azimuthal dependance 

As opposed to the axisymmetric energy kick imparted to 
our test-particle simulation, the merging satellite gala:xy is 
expected to impulsively heat the disc in a preferential direc- 
tion associated with its orbit. Therefore, we would expect 
to observe a phase difference in the distribution of substruc- 
ture observed in solar neighbourhood-like volumes located 
at different azimuthal angles. To explore this, we follow in 
our / = 20% simulation the E-L^ distribution of the parti- 
cles enclosed within four different spheres for ~ f Gyr. The 
spheres are located at (±8, 0, 0) kpc and (0, ±8, 0) kpc and 
hereafter we will refer to them as a; = ±8 kpc and y = ±8 
kpc. For this particular analysis we keep the sphere fixed 
at their initial locations. The results are shown in Figure [S] 
Each row of panels shows the time evolution of this dis- 
tribution on a different sphere. The panels are labelled, in 
the leftmost column, according to the position of the corre- 
sponding sphere. At any given time, the distribution of sub- 
structure in E-Lz space is different from sphere to sphere. 
A visual inspection of the first column of panels indicates a 
larger fraction of substructure on the y — —8 kpc sphere, 
0.1 Gyr after tmerger. Note that this direction corresponds 
with the direction of the first close passage of the satellite 
through the galactic disc, as shown by Figures 1 and 2 in 
Vff08. ft is also noticeable that the distribution of particles 
repeats itself along diagonals starting from bottom left to 
top right panels. This implies that the phase difference of 

the density wave can be well described as an spiral pattern, 

3 

travelling approximately an angle of -vr radians in 0.2 Gyr, 
in agreement with the spatial distribution of disc particles 
observed in Figure [3] To examine this, let us estimate the 
averaged angular period of a particle located at 8 kpc from 
the galactic centre. From Figure 2] at r = 8 kpc, (v) ~ 130 
km s~^. Thus, 



2nr 



: 0.38 Gyr, 



implying that in 0.2 Gyr, particles have travelled approxi- 
mately an angle tt through the disc. As a consequence, the 
ratio between the angular and the radial frequencie^ of the 
particles located inside our sphere is expected to be 

Figure [9] shows the distribution of orbital frequencies of 
the particles located inside the x — +8 kpc sphere, 0.7 Gyr 
af ter tmerger- The frequen cies were computed as described 
in I Gomez fc Helmil (|2010l . hereafter GHIO) and we refer the 
reader to their work for a detailed description of the ir esti- 
mation method (see also ICarpintero &: Aguilailll998l ). Note 
that particles tend to lie on the curve f20 = 2.1f2r-/3, in good 
agreement with our previous discussion. It is very interest- 
ing to observe that, in frequency space, particles are also 

Note that most of these particles have very low eccentricities 
and thus k Qr, where k represents the epyciclic frequency 
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Figure 5. Distribution of disc particles inside a sphere of 2.5 kpc 
radius located at 8 kpc from the centre, obtained from our A'-body 
simulation with a mass ratio of / = 20%, 0.7 Gyr after tmcrgcr. 
The left panel shows the distribution of particles in velocity space 
whereas the right panel corresponds to the E-L^ projection. Note 
that, unlike in velocity space, a large number of substructures can 
be observe in E-Lz space. 



distributed in well-defined lumps representing the different 
density waves crossing the volume. 



4.3 Dependance on satellite's mass 

It is well known that the amount of disc heating induced 
by a merger event strongly depends on the properties of the 
merging satellite, s uch as its conc entration and tot al mass 
(sec, e.g., Quin n. Hernguist fc FuUa gar 1993; Bcnso rTet al.l 



12004 iDi Matteo et all 12011 
1, mergers with 
to occur in 



VH08). Ahhough, since z = 
ratios as large as 20% are likely 
the likeli- 



mass 

galaxies such as the Milky Way, 
ness of t hese events inc reases with decreasing mass ratios 
(see e.g. iDe Lucia fc He lmi 2008; Kazantzidis et al. 20(3; 

to study 



ICooper et al.l bOloD . It is therefore interesting 



whether perturbations induced by less massive minor merg- 
ers can leave imprints on the phase-space distribution of 
galactic disc stars. In this section we investigate minor 
merger events with mass ratios of / = 15% and 10%. Ta- 
ble [T] summarises the main properties of both satellites. For 
less massive satellites than those considered here, dynami- 
cal friction is expected to be less efficient. Thus, this kind 
of satellites are unlikely to reach the centre of the host sys- 
tem and to cause significant changes in its structu re (VH08, 
iBinnev fc Tremaind[20oi : iKazantzidis et al.l \200^ . 

Figure [TD] shows, for both simulations, the time evo- 
lution of the E-Lz distribution of particles located inside 
a sphere centred at 8 kpc from the galactic central. Both 
spheres are rotating with an angular frequency derived from 
the corresponding velocity of their LSR. Interestingly, even 
for the simulation with a mass ratio / — 10% (bottom pan- 
els) , substructure can be clearly observed in the phase-space 
distribution of disc stars for a period of about 4 Gyr. This is 
approximately the same period of time that we were able 
to follow substructure in our / — 20% simulation. This 
can be explained as follows: Although the strength of the 
perturbation on the host's disc depends on satellite's mass, 
lower-mass satellites take longer to completely sink onto the 

* Note that, as shown in Figure 14 of VH08, the final thick disc's 
scale length of the / = 10% simulation does not significantly 
differs from the one measured on the / = 20% simulation. 
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galactic centre. For example, in our / — 10% simulation 
tmorgcr ~ 4 Gyr. As a consequence, a lower mass satellite 
can excite the galactic d isc for a longer period of time (see 
also iQuillen et a"l]|2009h . Note however that, after tmergcr, 
substructure can be followed only for ~ 2.3 and ~ 0.9 Gyr 
for the simulation with / = 15% and 10%, respectively. Low 
mass satellites generates density waves with smaller ampli- 
tudes and thus, their signatures after fmorgcr disappear faster 
as they quickly get within the noise level. 

4.4 Contamination from satellite's debris 

In the analysis carried out so far, only particles from the 
galactic disc were taken into account. As described by VH08, 
a large fraction of the particles from the satellite ends up in a 
disc-like configuration, with the same spatial orientation as 
the final thick disc. Therefore, the distribution of the satel- 
lite's particles in E-Lz space could partially overlap with 
the corresponding distribution of disc particles, complicat- 
ing the detection of density waves. 

The distribution of debris left by a satellite on a so- 
lar neighbourhood-like sphere, obtained from a simulation 
similar to the ones considered here, was already studied in 
detail by GHIO. The only difference between their / = 20% 
simulation and ours is the particle resolution of the satel- 
lite galaxy. To enhance the signatures of stellar streams, 
the satellite analysed by GHIO was modelled with a num- 
ber of particles ten times larger than the one used here (see 
also VH08). Figure 13 of GHIO shows the distribution of 
particles from both systems, in four different projections of 
phase-space. Not surprisingly, for a given energy, particles 
from the satellites tend to have smaller angular momentum 
than those from the disc (see bottom middle panel). Even 
in this simulation with a much higher satellite's resolution, 
substructure from the disc can be easily identified in E-Lz 
space. Note as well that substructure in the disc can be 
clearly observed in apocentre vs. pericentre space (bottom 
left panel). 



5 SUMMARY AND CONCLUSIONS 

In this work we have explored the impact that minor merg- 
ers have on the phase space distribution of a galactic stellar 
disc, such as the one in the Milky Way. Using numerical 
simulations, we showed that the pseudo-integrals of motion 
E-Lz space is well-suited to identify density waves in the 
disc's velocity field that have been generated by the merger 
event. In this space, particles in a given volume of physical 
space are distributed along features of nearly constant en- 
ergy, each associated with a different wave. The separation 
between density waves at different energy levels indicates the 
time since the disc was perturbed by the merging satellite. 
Interestingly, substructure in E-Lz space can be clearly ob- 
served even when solar neighborhood-like spheres with radii 
as large as 2.5 kpc are considered. This is in stark contrast 
with what is observed in the u-v plane, where the waves 
inside spheres with radii larger than ~ 0.25 kpc interfere, 
wiping out all signatures of ringing. An important implica- 
tion of this result is that a larger sample of stars could be 
explored to search for signatures of minor merger events in 
the Milky Way thick disc. This is relevant not only in the 




Figure 9. Distribution in frequency space of disc particles inside 
a sphere of 2.5 kpc radius located at 8 kpc from the galactic 
centre, obtained from our Af-body simulation with a mass ratio 
of / = 20%, 0.7 Gyr after fmorger- Note that substructure can be 
easily identify as well in this space. The solid line indicates the 
ratio between frequencies of disc particles, whereas the dashed 
lines denote the limits of this space. 

context of the forthcoming full 6D phase space Gaia cata- 
log, but also for available catalogs, suc h as those p rovided by 
RAVE (see, e.g., iBreddels et"ani2O10D and SDSS (|Lee et all 

(ml). 

By projecting the distribution of particles located inside 
spheres of 2.5 kpc radius onto E-Lz space we were able, for 
the first time, to study the properties of ringing in fully self- 
consistent A/^-body simulations of disc galaxies minor merg- 
ers. We showed that, in a simulation of a merger between a 
satellite galaxy and a host disc with a mass ratio / — 20%, 
signatures of ringing can be clearly observed for a total of 
~ 4 Gyr and for about 2.7 Gyr after the merger has con- 
cluded. It is remarkable that the < 1000 particles enclosed 
in our local volumes were enough to resolve a wealth of sub- 
structure in the resulting thick disc. In addition, we showed 
that the in-fall trajectory of the perturber satellite gives rise 
to an azimuthal dependence of the structure observed in E- 
Lz space. Note that we have only explored distributions of 
particles in this pseudo-integrals of motion space. Informa- 
tion about orbital properties of the merging satellite, such 
as its in-fall trajectory or its impact parameter, could be 
obtained by exploring the distributions of particles' angular 
phases as a function of galactocentric azimuthal angle. This 
will be analyse in a follow-up study. 

We have also explored the impact that mergers of dif- 
ferent total mass ratios have on the phase-space structure of 
the disc. Interestingly, mergers with mass ratios as small as 
/ — 10% give rise to density waves in the velocity field of the 
disc that can also be followed for a total of ~ 4 Gyr. This is a 
consequence of the longer period of time required by a lower 
mass satellite to fully merge with its host. On the other 
hand, once the merger concludes substructure can be fol- 
lowed for only ~ 0.9 and ~ 2.3 Gyr for simulations with mass 
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Figure 6. Time development of the E-Lz space for a distribution of particles located within a sphere of 2.5 kpc centred at 8 kpc from 
the galactic centre, obtained from our A^-body simulation with a mass ratio of / = 20%. The sphere is rotating with an angular velocity 
set by the velocity of the LSR. Note that, as expected, the number of lumps associated with density waves increases as time passes by. 
Substructure can be clearly observed for a total of ~ 4 Gyr. In this simulation the merger concludes at tmerger ~ 2 Gyr. 



-r 



^ -11 



t = 0.8 Gyr 

/ 


t= 1.1 Gyr 

/ 


t=1.5 Gyr 


t = 1 .8 Gyr 

/ 


t = 2.1 Gyr 

/ 


t = 2.4 Gyr 


t = 2.7 Gyr 


t = 3.0 Gyr 

•X- 


t = 3.3 Gyr 


t = 3.6 Gyr 
"^^^^ 


t = 3.9 Gyr 


t = 4.6 Gyr 



500 1500 

L^lkpc km/s] 



500 

Lj[kpc 



1500 

km / s] 



500 1 500 

L^lkpc km/s] 



500 1500 
L^lkpc km/s] 



500 

Lj[kpc 



1500 

km / s] 



500 1 500 

L^lkpc km/s] 



Figure 7. As in Figure [6] for particles located within a sphere centred at 5.47 kpc from the galactic centre. Note that mixing occurs 
more rapidly than in our sphere centred at 8 kpc, reflecting the dependance of the mixing time-scale on galactocentric radius. In this 



sphere, substructure can be clearly observe for a total of to ~ 2.8 Gyr. In this simulation the merger concludes at t,, 



: 2 Gyr. 



ratio of / = 10 and 15%, respectively. Note that these results 
are constrained by the low particle resolution of our A'^-body 
simulations; higher resolution simulations could allow us to 
follow substructure for longer periods of time. Thus, signa- 
tures of minor merger that have evolved for longer than 4 
Gyr could be identified from the very large and accurate 
sample of thick disc stars that soon will become available 
thanks to Gaia. 

M09 showed that a merger event ~ 1.9 Gyr ago could 
not only explain the high velocity streams seen in the 
Solar Neigh bour hood at approximate v of -60 (HR 1614, 
lEggenlll992l '). -80 jArifvanto fc Fuchsl200^). -1 00 (Arcturus), 
and -160 km/s (|Klement, Fuchs fc Ri^lHoOSl ') but also pre- 
dicts four new high velocity streams at v ~ 140, -120, 
40 and 60 km/s. By selecting samples of t hick disc stars 
(INordstrom et al ] |2004l : ISchuster' et. al ] |2006l ) they found in- 
dications of these new streams in velocity space. However, 
their results are not statistically significant since, due to 



the constraints on the size of the local volume (dmax = 
0.08 & 0.15 kpc), samples of only 451 and 766 stars were ex- 
plored. We are currently analysing the aforementioned 6D 
phase space catalogs from RAVE and SDSS to search for 
signatures of these streams in E-Lz space. 
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Figure 8. Azimuthal depcndance on the distribution of lumps in E-Lz space. Each row of panel shows the time development of the 
distribution of disc particles inside a different sphere, obtained from our A'^-body simulation with a mass ratio of / = 20%. The location 
of the spheres are indicated on the panel's bottom left corner of each panel. The variable 8 approximately indicates the angle with respect 
to the direction of the first close passage of the satellite through plane of the disc. The spheres are kept fixed in their initial locations. 
Note that, along diagonals starting from bottom left to top right panels, the distribution of particles approximately repeats itself. This 
indicates both, a phase difference on the density waves at different azimuthal angles and a relation between angular and radial frequencies 
of disc particles. 
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